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Summary 


A  Fortran  program,  written  by  A.G.P.  Warham,  for  calculating 
the  dynamics  and  acoustics  of  a  thin  plate,  has  been  extended  to 
include  excitations  other  than  the  unit  impulse.  Numerical  results 
of  acoustic  radiation  demonstrate  the  usefulness  of  the  program. 

In  the  acoustic  far-field  the  transient  sound  radiation  is  better 
computed  from  a  Fourier  integral  which  can  be  interpreted  in  terms 
of  drive  point  and  surface  radiation,  and  evaluated  numerically  by 
an  application  of  the  fast  Fourier  transform  algorithm.  Numerical 
examples  are  given  using  both  'thin'  and  ' thick; •  plate  theories. 
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1.  INTRODUCTION 


The  infinite  fluid-loaded  plate  is  one  of  the  prototype 
problems  in  the  field  of  fluid-structure  interaction,  and  its 
time-harmonic  response  and  associated  sound  radiation  have  been 
studied  extensively  in  the  acoustic  literature.  The  standard 
texbook  of  Junger  &  Feit  C1I  gives  the  mathematics  and  the 
physics,  and  Spicer's  C2D  intensity  vector  plots  are  vivid 
illustrations  of  the  main  features  of  the  acoustic  field.  The 
transient  dynamics  and  acoustics  of  the  infinite  fluid-loaded 
plate  have  not  received  nearly  so  much  attention  despite  their 
importance  as  canonical  problems  in  underwater  acoustics. 

Warham  C33  has  recently  written  a  Fortran  computer  program 
for  calculating  the  response  of,  and  sound  radiation  from,  an 
infinite  thin  plate  which  is  excited  by  a  unit  point  impulse. 

The  program  cleverly  combines  analytical  and  numerical  approxi¬ 
mations  to  achieve  efficient  and  accurate  computations  over  a 
wide  range  of  geometric  values  and  material  properties.  Many 
references  to  previous  theoretical  work  can  be  found  in  Warham's 
comprehensive  report. 

In  the  small  amount  of  time  set  aside  for  this  work  it  was 
the  aim  to  (a)  extend  the  range  of  applicability  of  Warham's 
program  by  giving  the  user  a  choice  of  transient  force  excita¬ 
tions;  (b)  give  some  numerical  results  of  sound  radiation  which 
demonstrate  that  Warham's  program  is  a  useful  contribution  to 
research  in  transient  fluid-structure  interaction;  (c)  test  a 
stationary  phase  approximation  to  the  transient  pressure  in  the 
far-field;  and  (d)  take  a  first  stumble  in  the  direction  of 
developing  a  software  capability  for  making  predictions  of 
transient  sound  radiation  from,  and  transient  sound  scattering 
by,  submerged  elastic  structures. 


2.  PROBLEM  FORMULATION 


An  infinite  plate,  see  Figure  1,  lies  in  the  x,y  plane. 

The  upper  half space  contains  a  fluid  and  the  lower  half space  a 
vacuum.  The  plate  is  excited  by  a  transient  stress  distribution 
F(r,t),  the  plate's  normal  displacement  W(r,t)  in  bending  and 
the  acoustic  pressure  p(r,z,t)=p(R,0,t)  being  required. 

The  plate's  displacement  is  assumed  to  satisfy  the  thin 
plate  theory  of  Sophie  Germain,  viz. 

DCa2/ar2-*-(l/r)3/3rD2W(r,t)+M32W(r,t)/at2  *  F(r,t)-p(r,0,t> 

(2.1) 

where  D=Eh2/12( l-o2)  is  the  plate's  flexural  rigidity  in  which 

E  is  the  Young's  modulus,  o  is  the  Poisson's  ratio  and  h  is  the 

thickness;  M=p  h  is  the  plate's  mass  per  unit  area  in  which 

s  s 

is  the  density.  The  acoustic  pressure  p(r,z,t)  satisfies  the 

wave  equation 

Ca2/ar2+( l/r)a/3r+a2/az23p(r,z,t)  *  a2p(r,z,t)/c2at2  (2.2) 
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in  which  c  is  the  velocity  of  sound  in  the  fluid. 

The  solution  to  the  above  equations  proceeds  by  defining 
the  Fourier-Bessel  transform  pair  for  a  quantity  U(r,z,t),  i.e. 

U(r,z,t)=  (l/2ir)/^  /*  tJ(a,z,w)J0(ar)exp(-iujt)o<dadw 

(2.3) 

U(a,z,u>)  =  /"  U(r,z,t)J0(<xr)exp(iut)rdrdt 

Thus,  after  defining  the  plate's  displacement  W(r,t)  and  the 
pressure  p(r,z,t)  by  the  transforms 

W(r,t)=  (l/2ir)/^  /*  W(a,w) JQ(ar )exp(-iu)t)odotdui  (2.4) 

p(r,z,t)=  (l/2ir)/^  /“  p(a,z,u) JQ(ar )exp( -iwt)otdadij)  (2.5) 

it  is  possible  to  show  that 

W(a,u>)  =  F(cx,w)  /  (Da*-uj2M-ipu)2/Y)  (2.6) 

p(a,z,ijj)  *  -ipw2-y  1F(a,w)exp(  iyz )  / (Da4-iu2M-ipw2/Y)  (2.7) 

in  which  F(a,u>)  is  the  transform  of  the  transient  stress  exci- 
2  2  1/2 

tation;  y=+(k  -a  )  in  order  to  satisfy  the  radiation  condition 
p  is  the  fluid's  density;  and  k  is  the  acoustic  wavenumber,  w/c. 

The  acoustic  momentum  relation  3p(r,z,t ) /3z=-p32W(r ,z,t) /3t2  has 
been  used  to  ensure  continuity  of  the  plate's  displacement  W(r,t) 
and  the  acoustic  particle  displacement  W(r,z,t)  at  the  plate's 
surface,  z=0. 

The  integrals  in  equations  (2. 4 -2. 5)  must  be  evaluated 
numerically  because  closed-form  expressions  are  not  available. 
Mackertich  &  Hayek  C43  have  given  analytical  approximations  to 
the  acoustic  pressure  for  the  specific  case  of  impulse  exci¬ 
tation;  however,  according  to  Warham  some  of  their  approximations 
are  lacking  in  accuracy.  Despite  this  (alleged)  shortcoming, 
Makertich  and  Hayek  have  identified  correctly  much  of  the  physics 
of  this  transient  sound  radiation  problem. 
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an  excitation  which  is  independent  of  the  wavenumber  a. 

(b)  Impulse 

The  impulse  applied  at  t  =  0  is  defined  as 

Q(t)  -  P6(t>  (3.3) 

where  P  is  the  magnitude  of  the  impulse.  Its  transform  is 
obtained  from  equation  (3.2)  as 

F(u>)  =  P/2w  (3.4) 

This  is  the  excitation  used  in  Warham's  original  Fortran  program, 
which  has  now  been  extended  to  include  the  transient  excitations 
described  in  (c)-(g)  below. 

( c )  Sine  Wave 

The  sine  wave  which  is  switched  on  at  t=0  and  switched  off 
at  t=T  is  defined  as 

Q(t)  =  FQsin(u»0t),  for  0<t<T  (3.5) 

Q(t)  =  0,  for  t<0  and  t>T 

the  time  for  one  cycle  being  l/fQ  =  2ir/u0.  Its  transform  is 

F(w)  =  (F0/2x)Cexp(  iwT)  £iujsin(w0T) -w0cos(u0T)  1+WqD/ (Wq-oj2  ) 

(3.6) 

This  excitation  taken  over  a  half -cycle  is  sometimes  used  to 
simulate  the  force  due  to  an  elastic  impactor. 

(d)  Squared  Sine  Wave 

The  squared  sine  wave  which  is  switched  on  at  t=0  and  off 
at  t=T  is  defined  as 

Q(t)  =  F0sin2(wQt),  for  0<t<T  (3.7) 

Q(t)  =  0,  for  t<0  and  t>T 

the  time  for  a  half-cycle  (one  loop)  being  l/2f0=w/w0.  Its 
transform  is 

F(u>)  =  <F0/4*iw)Cexp<  iwT)-13  (3.8) 

- ( FQ/4ir )  Cexp(  iwT)  Ciwcos  (u)QT)  +2w0sin(  2u»QT)  }-iw3/  ( 4w2-w2 ) 

This  excitation,  taken  over  a  half -cycle,  has  been  recommended  by 
Akay  &  Latcha  C5D  for  simulating  an  elastic  impact,  because  it 
gives  a  smoother  acceleration  response  than  does  the  sine  wave. 


(e)  Attenuated  Sine  Wave 


The  attenuated  sine  wave  which  is  switched  on  at  t=0  and  off 
at  t=T  is  defined  as 

Q(t)  =  F0exp(-at)sin(u)0t) ,  for  0<t<T  (3.9) 

Q(t)  =  0,  for  t<0  and  t>T 

the  time  for  one  cycle  being  1/f 0=2tt/uj0.  Its  transform  is 

F(w)  =  (FQ/2ir)Cexp(ixT)CiXsin(w0T)-u)0cos(u0T>  3+WqD/ (Wq-X2) 

(3.10) 

in  which  X  =  w+ia.  This  is  the  excitation  used  by  Bolgov  & 
Nikiforov  C63  in  their  analysis  of  the  impact  excited  infinite 
plate. 

( f )  Triangular  Pulse 

The  triangular  pulse  which  is  switched  on  at  t=0  and  off 
at  t=T  is  defined  as 

Q(t)  =  ( 2FQ/T) t ,  for  t=0<t<T/2 

Q < t )  =  - ( 2FQ/T) t+2FQ ,  for  T/2<t<T  (3.11) 

Q(t)  =  0,  for  t<0  and  t>T 

Its  transform  is 

F(w)  =  (-2F0/2Tru)2T)Cl-exp(iwT/2)32  (3.12) 

Nakayama  &  Nakamura  C7]  have  used  this  pulse  as  an  acoustic 
excitation  in  their  work  on  transient  sound  radiation  from 
circular  plates. 

(g)  Rectangular  Pulse 

The  square  pulse  which  is  switched  on  at  t=0  and  off  at 
t=T  is  defined  as 

Q(t)  =  Fq,  for  0<t<T  (3.13) 

Q(t)  *  0,  for  t<0  and  t>T 

Its  transform  is 

F(u)  *  (FQ/2iriu;)Cexp(iwT)-13  (3.14) 

This  pulse  can  be  used  to  simulate  a  finite  width  Impulse.  It 
is  also  the  basic  component  of  extended  rectangular  waveforms, 
such  as  the  finite  length  Walsh  functions  and  pseudo  random 
binary  sequences.  It  would  not  be  difficult  to  add  these  two 
waveforms  to  the  list,  as  their  Fourier  transforms  are  known. 
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4. 


APPROXIMATION  TO  FAR -FIELD  PRESSURE 


(a)  Stationary  Phase 

If  the  a-integration  in  equation  (2.5)  is  done  by  the 
method  of  stationary  phase  C13,  then  the  the  far-field  pressure 


pf  <R,0,t)  =  (-p/2irR) uj2CF(<x0 ,w)  /G(aQ ,w) 3expC-iw(  t-R/c ) 3du> 


(4.1) 


results,  where 


G(<Xq,u)  =  Da4-w2M-ipu)2/YQ 

aQ=k.sin(0)  (4.2) 

Y0=k. cos(8) 


The  integral,  which  is  invariant  with  respect  to  the  factor 
t-R/c,  can  be  evaluated  simultaneously  for  equally  spaced  time 
intervals  by  an  application  of  the  fast  Fourier  transform  algor¬ 
ithm.  Its  use  for  transient  sound  radiation  calculations  does 
not  appear  to  have  been  explored  in  any  detail,  in  the  acoustics 
literature. 

(b)  Physical  Interpretation 

The  integration  in  equation  (4.1)  can  be  done  by  Cauchy's 
residue  theorem,  as  there  are  no  branch  points  to  'muddy  the 
waters'.  The  function  F(a0,u)  will  introduce  some  complications 

into  the  analysis,  but  it  is  possible  to  proceed  in  a  cavalier 
way  by  considering  only  the  contribution  made  by  the  function 
G(a0,uj),  which  is  responsible  for  simple  poles  whose  locations 

are  the  roots  of  the  reduced  cubic  equation 

a0+ajW+a3uj3=  0  (4.3) 

where 

aQ  =  -ipc/cos0 

ax  =  -M  (4.4) 

a3  =  (D/c4)sin4© 

At  an  arbitrary  angle  0  there  are,  in  the  absence  of 
structural  damping,  three  roots  of  the  farm 

*  -2iw^ 

w2  *  wr+*wi 

u>3  *  -u»r+iwi  (4.5) 

in  which  u>r  and  are  positive. 

For  a  5cm  thick  steel  plate  with  water  above  and  a  vacuum 
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below,  with  constants  as  defined  in  Section  5  below,  the  values 
of  u>r  and  u>i  are 


e 

wrxl0  3 

wixl0-3 

fr 

10 

983 

2.0 

150.1 

20 

253 

2.1 

40.3 

30 

119 

2.2 

18.9 

40 

72 

2.5 

11.4 

50 

51 

3.0 

8.1 

60 

40 

3.8 

6.4 

70 

35 

5.2 

5.5 

00 

34 

8.5 

5.4 

85 

37 

12.8 

5.9 

89 

55 

26.5 

8.7 

where  fr=u>r/2ir  is  the  oscillating  frequency,  in  kHz.  It  is  easy 
to  show  that  the  roots  may  be  approximated  by  the  formulae 

u>r  =  c2(M/D)1/2/sin20 

uk  =  pc/(2Mcos0)  (4.6) 

for  values  of  0  less  than  70°,  which  shows  that  while  u)r  depends 

on  both  the  mass  and  stiffness  properties  of  the  plate,  uk  is 

controlled  by  its  mass  properties.  Enthusiasts  of  infinite  plate 
problems  will  recognise  the  minimum  value  of  in  the  approxi- 

2  1/2 

mation,  viz.  wc=c  (M/D)  ,  as  the  coincidence  frequency,  which 

can  be  defined  as  the  frequency  at  which  the  phase  velocity  of 
straight-crested  waves,  on  an  unloaded  plate,  has  the  value  of  the 
speed  of  sound  in  the  fluid.  This  frequency  is  4.7kHz  for  a  5cm 
thick  steel  plate  in  water.  At  frequencies  higher  than  the  coinci¬ 
dence  frequency  there  is  a  strong  lobe  in  the  radiated  field  when 
time-harmonic  conditions  prevail,  the  angle  at  which  the  lobe 

-1  1/2 

occurs  being  about  6c  =  sin  (u>c/w)  .  Thus,  it  is  possible  to 

assign  meaning  to  as  the  true  frequency  of  coincidence  lobe 

radiation  at  the  angle  0.  See  Junger  &  Feit  C13  for  a  discussion 
of  coincidence  lobe  radiation. 

Three  contours  of  integration  are  necessary  in  order  to 
evaluate  the  Fourier  integral  analytically: 

(1)  The  first,  for  t<0,  is  the  semicircle  in  the  upper 
half  plane,  the  diameter  of  which  runs  along  the  real  u> 
axis,  indenting  the  poles  u>2  and  u>3.  Thus,  because 

there  are  no  poles  enclosed  by  this  contour  the 
pressure  is  zero  for  t<0,  which  is  consistent  with 
causality. 
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(2)  The  second,  for  t=0  to  R/c,  is  the  semicircle  in 
the  upper  half plane,  the  diameter  of  which  runs  along 
the  real  w  axis.  Thus  the  poles  t>  2  and  w3  are  enclosed. 

by  the  contour.  These  poles  have  no  separate  existence; 
rather,  they  combine  to  give  a  real  value  of  pressure 
which  is  proportional  to  expCw^ t-R/c )3sinCwr( t-R/c )+#3. 

Bearing  in  mind  that  R  is  infinite  in  the  far-field, 
this  is  a  pressure  whose  amplitude  increases  exponen¬ 
tially  with  time,  from  zero  to  a  maximum  at  t=R/c,  and 
which  oscillates  at  a  radian  frequency  of  w  . 


(3)  The  third,  for  t>R/c,  is  the  semicircle  in  the 
lower  half plane,  the  diameter  of  which  runs  along  the 
real  w  axis.  Thus  only  the  pole  is  enclosed,  giv¬ 
ing  an  exponentially  decreasing  pressure  proportional 
to  expC-2wi ( t-R/c ) 3 . 


When  9=0°  (normal  to  plate)  the  coefficient  a3  vanishes, 
and  the  only  root  is 

=  -i(pc/M)  (4.7) 

which  lies  in  the  lower  halfplane  and  gives  an  exponentially 
decaying  pressure  field  proportional  to  expC-( pc/M) ( t-R/c ) 3 ,  for 
t>R/c.  Because  there  are  no  roots  in  the  upper  halfplane  the 
pressure  is  zero  for  t<R/c.  This  root  is  independent  of  the 
plate's  stiffness,  D,  depending  only  on  the  fluid's  constants 
and  the  plate's  mass  per  unit  area. 

The  aforementioned  behaviour  of  the  roots  can  be  interpreted 
physically  with  the  help  of  the  diagram 


d 


in  which  it  is  posible  to  show  that  R^=R-d.sin0  in  the  acoustic 

far-field.  Now  consider  an  energy  packet  that  is  injected  into 
the  plate  at  time  t=0:  some  of  it  will  propagate  to  P  via  the 
direct  path  (OP)  in  time  R/c,  and  some  of  it  will  propagate  to  P 

via  the  indirect  path  (OQP)  in  time  d/c  +R/c-d. sin0/c,  where  c 

g  g 

is  the  group  velocity  of  the  wave  packet,  in  the  plate,  which  is 
losing  energy  by  radiation.  When  these  two  travel  times  are 
subtracted  it  becomes  evident  that  it  is  possible  for  energy  to 
arrive  at  P  by  the  indirect  path  OQP  before  it  arrives  by  the 
direct  path  OP,  provided  c^>c/sin0.  This  is  the  condition, 

assuming  the  group  velocity  is  greater  or  equal  to  the  phase 
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velocity,  for  radiation,  at  an  angle  9,  by  supersonic  straight- 
crested  waves  travelling  on  a  fluid-loaded  plate.  Thus,  super¬ 
sonic  waves,  on  the  plate,  are  the  mechanism  by  which  radiation 
prior  to  time  t=R/c  is  produced.  See  Junger  &  Feit  for  general 
details  of  radiation  from  wave-bearing  surfaces,  and  Crighton  C8D 
for  a  full  discussion  of  waves  on  a  fluid-loaded  plate. 

This  non-rigorous  treatment  of  the  Fourier  integral  has 
shown  that  the  roots  of  the  function  G(ot,Q,u))  have  a  physical 

significance;  the  first  root  can  be  identified  with  radiation  from 
the  excitation  point,  because  it  only  contributes  from  time  t=R/c, 
the  direct  arrival  time;  the  other  two  roots  (in  combination) 
can  be  identified  with  surface  radiation  because  they  only  con¬ 
tribute  before  time  t=R/c,  which  is  only  possible  if  energy  has 
travelled  supersonically  along  the  plate  surface  before  radiating 
into  the  fluid.  Now,  as  the  distance  d  (see  diagram  above)  along 
the  plate  increases,  indirect  radiation  from  the  plate  arrives 
increasingly  earlier  than  the  direct  radiation.  Furthermore,  as 
d  increases,  the  amplitude  of  the  supersonic  wave  in  the  plate 
will  have  decreased  exponentially  because  of  radiation  and  struc¬ 
tural  damping,  which  means  that  its  radiation  will  also  decrease 
exponentially,  assuming  that  it  is  proportional  to  the  the  wave's 
amplitude.  As  the  distance  d  reduces,  it  follows  that  the 
indirect  radiation  will  increase  exponentially  until  the  origin 
is  reached,  at  which  point  it  merges  continuously  with  the  direct 
radiation. 

Precise  details  of  the  far-field  pressure  pulse  will  depend 
on  the  frequency  spectrum  of  the  excitation,  but  general  obser¬ 
vations  are  possible  for  the  case  of  an  excitation  switched  on  at 
time  t=0  and  then  switched  off  at  time  t=T.  Before  the  time  of 
the  direct  arrival,  t=R/c,  the  pressure  will  oscillate  at  a 
frequency  whose  value  is  dependent  on  the  angle  of  observation, 
the  minimum  value  being  close  to  the  coincidence  frequency.  The 
amplitude  of  this  pressure  will  increase  exponentially  with  time 
until  t=R/c,  at  which  time  it  will  merge  continuously  with  the 
direct  pressure.  For  times  t=R/c  to  t=T+R/c  the  pressure  will  be 
the  sum  of  a  relatively  smooth  direct  wave  whose  travel  time  is 
t=R/c,  and  the  oscillating  indirect  wave  which  will  appear  at 
times  to  anticipate  sharp  discontinuities  in  the  excitation.  For 
times  t>T+R/c  the  pressure  field  will  smoothly  decay  to  zero. 

When  9=0°  the  pressure  is  zero  before  t=R/c;  it  is  due  only  to 
direct  radiation  at  times  t=R/c  to  T+R/c ;  and  it  decays  smoothly 
to  zero  for  t>T+R/c.  This  description  of  the  radiated  pulse  does 
not  permit  the  'ringing',  after  time  t=T+R/c,  which  was  predicted 
by  Mackertich  &  Hayek;  thus  it  supports  Warham's  view  that  their 
results  are  in  error. 

Warham  gives  explicit  formulae  for  the  residue  contri¬ 
butions,  for  the  case  of  impulse  excitation,  and  there  is  agree¬ 
ment  in  those  areas  of  the  analysis  where  overlap  occurs.  It 
is  to  be  hoped  that  a  proper  treatment  of  the  stationary  phase 
approximation  and  the  Fourier  integral  evaluation  by  contour 
integration  will  be  undertaken  and  reported  later. 

(c)  Timoshenko-Mindlin  Plate  Theory 

An  advantage  of  using  the  fast  Fourier  transform  algorithm 
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to  evaluate  the  integral  in  equation  (4.1)  is  that  different 
plate  theories  may  be  used  without  the  complication  of  finding 
the  roots  of  the  function  G(otQ,u>),  which  is  a  simple  polynomial 

only  for  Sophie  Germain  plate  theory. 

For  Timoshenko-Mindlin  plate  theory 

G(<x0,u>)  =  CC(DaQ-psh3w2/12)  (a^-psw2/KG)-psha)2}/ 

(l+Da2/KGh-pgh2w2/12KG)3  -  ipw2/y0  (4.8) 

in  which  G=E/2(l+a)  and  k=10 ( 1+a) / ( 12+lla) .  For  exact  linear 
plate  theory  the  function  G(<x0,u)  can  be  found  from  the  work  of 

Spicer  C93. 

5.  NUMERICAL  RESULTS 

(a)  Constants 

The  numerical  results  of  the  transient  acoustics  of  a 
steel  plate  in  water,  shown  in  Figures  2-11,  were  obtained  using 
the  following  material  and  geometric  constants: 

Steel  Plate:  Young's  modulus  19.5x10*° 

Poisson's  ratio  0.29 
Density  7700.0 
Thickness  0.05 
Hysteretic  loss-factor  0.01 

Water:  Density  1000.0 

(one  side  only)  Sound  velocity  1500.0 

The  use  of  a  constant  hysteretic  loss-factor  makes  the  responses 
non-causal  (see  Warham  for  a  discussion),  but  this  is  barely 
noticeable  in  Figures  2-11  .  On  each  plot,  where  applicable,  are 
marked  the  times  of  the  direct  travel  time  of  the  leading  edge  of 
the  excitation,  D=R/c,  the  minimum  indirect  acoustic  travel  time, 
I=R.cos8/c,  and  the  direct  travel  time  for  the  trailing  edge  of 
the  excitation,  L=T+R/c. 

(b)  Near-field  Acoustic  Pressure 

Figures  2-5  show  transient  sound  radiation  versus  time, 
obtained  from  Warham' s  program,  for  various  excitations.  In 
each  plot  the  angle  of  measurement  is  8=45°.  Because  the 
pressures  are  in  the  'near-field',  the  physical  interpretations 
of  Section  4  will  be  incomplete  since  the  domains  of  the  subson¬ 
ic  wave  and  the  supersonic  'leaky'  wave  have  been  entered  -  see 
Crighton  C83.  Nonetheless,  the  results  are  expected  to  show 
certain  features  that  are  present  in  the  far-field,  viz.  an 
oscillating  indirect  arrival  before  t=R/c,  and  a  smooth  relax¬ 
ation  to  zero  after  the  time  of  arrival  t=T+R/C  of  the  trailing 
edge  of  the  pulse.  In  the  far-field  the  indirect  arrival  oscil¬ 
lates  at  the  frequency  wr,  which  has  been  defined  as  the  true 

coincidence  frequency  at  the  angle  of  measurement  9,  and  it  has 
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no  well  defined  beginning.  In  the  near-field  the  nature  of  the 
indirect  arrival  can  be  surmised  with  the  help  of  the  diagram  of 
Section  4b,  in  which  the  point  P  is  now  in  the  near-field.  The 
first  acoustic  arrival  cannot  occur  earlier  than  time  t=R.cos0/c, 
via  the  path  OQP  in  which  P  is  directly  above  Q;  this  arrival  is 
associated  with  an  infinite  (at  least  mathematically)  value  of 
wr.  The  last  indirect  acoustic  arrival  is  at  time  t=R/c,  via 

the  path  OQP  in  which  0  and  Q  are  coincident;  this  arrival  is 
associated  with  the  value  of  wr  that  corresponds  to  the  angle 

of  observation,  e.g.  wr=  9.5kHz  at  6=45°.  Thus,  the  indirect 

arrival  has  a  well  defined  beginning,  at  time  t=R.cos0/c  where 
the  oscillation  frequency  is  infinite,  and  ends  at  time  t=R/c 
with  an  oscillation  frequency  of  wr. 

In  Figure  2a  the  excitation  is  a  unit  impulse  whose 
frequency  spectrum  is  a  constant.  The  distance  R  is  1.0m.  From 
the  time  of  the  indirect  arrival,  the  pressure  increases  rapidly 
with  increasing  period  of  oscillation;  the  oscillations  are  too 
rapid  to  show  on  the  figure.  The  pressure  then  decreases  with 
increasing  period  of  oscillation  until  the  time  of  direct  travel 
when  it  almost  vanishes.  These  results  were  obtained  with  an 
upper  limit  of  frequency  integration  set  to  1MHz,  Warham's  pro¬ 
gram  providing  an  additional  contribution  by  means  of  an  analyti¬ 
cal  formula.  Figure  2b  shows  the  results  obtained  without  the 
extra  analytical  contribution.  The  non-causal  behaviour  is 
marked,  and  the  early  time  arrival  parts  of  the  plots  differ 
considerably.  When  the  upper  limit  of  integration  is  truncated 
at  10kHz,  the  maximum  level  (not  shown)  is  less  than  unity, 
which  implies  the  dominance  of  high  frequency  surface  radiation. 

As  a  test  to  see  whether  or  not  steady-state  conditions  are 
reached  quickly,  a  1kHz  and  a  10kHz  sine  wave  excitation  were 
switched  on  for  5  cycles  each.  The  results  are  shown  in  Figures 
3a  and  3b,  in  which  the  distance  R  is  1.0m.  The  plots  show  that 
after  an  Initial  transient,  time-harmonic  conditions  prevail 
within  1  cycle.  Spicer  C10J  has  written  a  Fortran  program  for 
calculating  near-field  pressures  due  to  time-harmonic  excitation. 
This  program  was  used  to  confirm  that  the  steady-state  levels 
shown  in  Figures  3a  and  3b  are  exactly  as  expected. 

In  order  to  illustrate  the  change  in  pressure  shape  as  a 
disturbance  propagates  from  the  near-field  to  the  far-field, 
a  square  wave  was  switched  on  for  1/2  millisecond.  Figure  4a 
shows  the  pressure  at  R=0.2m.  After  a  small  rapidly  changing 
initial  transient,  at  the  time  of  the  indirect  arrival,  the 
pressure  rises  rapidly  to  reach  a  maximum  close  to  the  time  of 
the  direct  arrival.  The  pressure  then  drops  smoothly  until  the 
indirect  arrival  part  of  the  'switch-off'  occurs  which  causes 
another  rapidly  varying  transient,  ending  when  the  pressure  drops 
sharply  to  a  negative  value  after  which  it  rises  smoothly  to  zero. 
It  is  believed  that  the  prevalent  feature  in  the  pressure  is  due 
to  acoustic  radiation  from  the  excitation  point,  and  the  rapidly 
varying  oscillations  are  due  to  radiation  from  the  plate  surface. 
This  view  is  reinforced  by  the  plots  shown  in  Figures  4b,  5a  and 
5b,  which  are  for  R*l,  5  and  25m  respectively.  The  shape  of  the 
pressure  from  the  drive  point  is  (possibly)  unaltered,  but  it 
is  masked  by  oscillations,  due  to  the  indirect  plate  radiation. 
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which  decrease  in  frequency  to  about  10kHz,  as  the  distance 
increases  to  25m.  According  to  the  Table  in  Section  4b,  the 
oscillations  will  have  decreased  to  9.5kHz  by  the  time  the 
acoustic  far-field  has  been  reached.  Figure  4b  is  a  vivid  illus¬ 
tration  of  the  pressure  apparently  anticipating  the  switching 
on  and  off  of  the  excitation,  which  feature  is  due,  of  course,  to 
indirect  plate  radiation. 

( c)  Far-Field  Sound  Radiation 


A  Fortran  program  has  been  written  to  calculate  the  far- 
field  radiation  by  evaluating  the  Fourier  integral  in  equation 
(4.1)  using  an  application  of  the  fast  Fourier  transform  algor¬ 
ithm.  Details  of  the  integration  scheme  chosen  are  not  given 
herein.  Both  Sophie  Germain  'thin'  and  Timoshenko-Mindlin 
'thick'  plate  theories  have  been  included  as  options  in  the 
computer  program.  Numerical  results  are  shown  in  Figure  6-11. 

The  levels  have  been  adjusted  to  the  reference  distance  of  lm, 
and  although  the  time  scale  is  arbitrary  it  is  useful  to  know 
that  the  the  direct  arrival  time  has  been  set  to  time  t=lms,  by 
the  program,  by  use  of  a  variable  RP  whose  value  is  in  the  graph 
caption. 

Figure  6  shows  transient  pressures  in  the  far-field  when 
the  excitation  is  a  unit  Impulse.  In  Figure  6a  the  angle  of 
observation  is  0=0°;  the  pressure  is  zero  for  time  t<R/c,  it 
is  basically  a  sharp  spike  at  time  t=R/c,  and  then  it  vanishes 
almost  immediately.  In  Figure  6b  the  angle  of  observation  is 
0=45°;  the  pressure  rises  exponentially,  with  a  9.5kHz  modulation, 
until  t=R/c,  at  which  time  it  relaxes  smoothly  to  zero.  These 
plots  are  in  agreement  with  the  physics  described  and  the  values 
tabulated  in  Section  4b. 

As  a  test  to  see  whether  or  not  steady-state  conditions 
are  reached  as  quickly  in  the  far-field  as  they  are  in  the  near¬ 
field,  1kHz  and  10kHz  sine  waves  were  switched  on  for  5  cycles 
each.  Figure  7a  shows  that  after  an  initial  transient,  oscil¬ 
lating  at  the  frequency  of  9.5kHz,  as  expected,  time-harmonic 
conditions  are  obtained  within  1  cycle  for  the  1kHz  excitation. 

In  Figure  7b  the  indirect  arrival  is  large  because  the  excitation 
frequency  is  close  to  the  frequency  of  its  oscillation.  Steady 
state  conditions  are  in  fact  reached  within  1  cycle,  but  this  is 
masked  by  the  indirect  arrival  of  the  transient  generated  by  the 
switch-off.  Calculations  using  standard  formulae,  found  in  Junger 
&  Feit,  have  confirmed  that  the  correct  steady-state  levels  have 
in  fact  been  achieved. 

In  Figures  8-9  the  excitation  is  a  rectangular  pulse  that 
has  been  switched  on  for  1/2  millisecond.  Figure  8a  shows  the 
pressure  at  normal  incidence,  0=0°:  the  pressure,  which  is  zero 
before  time  t=R/c,  rises  rapidly  as  the  leading  edge  reaches  the 
measuring  position,  drops  slowly  until  the  time  t=T+R/c  of  the 
switch  off  when  it  drops  rapidly  to  a  negative  value,  after  which 
it  rises  smoothly  to  zero.  In  Figure  8b  the  angle  of  measurement 
is  0=20°:  before  t=R/c  the  indirect  arrival  due  to  the  leading 
edge  of  the  excitation  is  responsible  for  the  exponentially 
increasing  pressure,  modulated  by  an  oscillation  of  about  40kHz; 
between  times  t=R/c  and  T+R/c  the  direct  arrival  is  modulated  by 
the  oscillations  of  the  indirect  arrival  due  to  the  trailing  edge 
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of  the  excitation;  for  times  t>T+R/c  the  pressure  rises  smoothly 
to  zero.  Figures  9a  and  9b,  for  9=45°  and  0=70°,  respectively, 
show  similar  features,  except  that  the  oscillation  rates  are 
about  9.5  and  5.5kHz.  The  physics  is  well  explained  by  the  text 
and  the  Table  in  Section  4.  A  comparison  of  Figures  5b  (Harham 
at  25m)  and  9a  (James  in  the  far-field)  shows  that  these  plots 
are  very  similar  in  shape,  and  in  level  if  Warham's  results  are 
referred  to  lm  using  a  spherical  spreading  correction,  which 
reinforces  confidence  in  the  correctness  of  the  Fortran  programs. 

Figure  10  and  11  are  comparisons  of  the  transient  pressure 
pulses  obtained  from  the  'thin'  and  'thick'  plate  theories  at 
0=45°.  In  Figure  10  the  excitation  is  a  rectangular  pulse  that 
has  been  switched  on  for  1/2  millisecond,  the  pressure  level 
obtained  by  using  'thick'  plate  theory  being  shown  in  Figure  10a 
and  the  pressure  level  obtained  by  using  'thin'  plate  theory  being 
shown  in  Figure  10b.  The  two  plots  are  qualitatively  similar;  but 
the  'thick'  plate  pressure  oscillates  at  a  higher  frequency  than 
the  'thin'  plate  pressure,  and  it  has  two  sharp  spikes  at  its 
leading  and  trailing  edges.  The  reason  for  the  former  feature  is 
that  the  frequency  ur  associated  with  a  given  aspect  0  is  higher 

for  'thick'  plate  than  it  is  for  'thin'  plate  theory.  The  reason 
for  the  latter  is  not  understood,  but  jump  discontinuities  in  the 
excitation  appear  to  produce  more  of  an  overshoot  in  the  acoustic 
pulse  when  'thick'  plate  theory  is  used.  Also,  the  Fourier  inte¬ 
gral  is  only  conditionally  convergent,  at  times  t=R/c  and  t=T+R/c, 
when  'thick'  plate  theory  is  used  and  the  excitation  spectrum 
drops  only  as  fast  as  1/w;  which  means  that  it  is  not  conclusive 
that  the  correct  answers  are  obtained  at  these  times.  In  Figure 
11  the  excitation  is  the  triangular  pulse  on  for  1/2  millisecond: 
the  plots  are  very  similar,  sharp  spikes  being  absent  when  'thick' 
plate  theory  is  used  because  there  are  no  jump  discontinuities  in 
the  excitation. 


6.  CONCLUDING  REMARKS 


The  limited  number  of  numerical  results  shown  have  dem¬ 
onstrated  that  Warham's  program  is  a  valuable  contribution  to 
research  in  transient  sound  radiation.  The  program  is  particu¬ 
larly  flexible  because  it  copes  with  the  acoustics  (and  the 
dynamics)  both  at  short  and  long  ranges.  Its  major  drawback  is 
its  size  of  almost  9000  lines,  and  its  reliance  on  library  (NAG) 
subroutines  that  may  not  be  readily  available.  In  addition  to 
providing  bench-mark  computations,  this  Fortran  program  could 
be  used  to  investigate: 

(a)  The  build-up  and  eventual  decay  of  an  acoustic 
energy  vortex  as  a  sine  wave  excitation  is  switched 
on,  for  a  time  sufficient  for  time-harmonic  response 
to  prevail,  and  then  switched  off.  Spicer's  C23 
plots  show  energy  vortices  and  Waterhouse  et  al  C113 
give  a  condition  for  their  occurrence. 

(b)  The  linear  dynamics  and  acoustics  due  to  impact 
by  an  elastic  sphere.  For  this  problem,  the  papers 
of  Akay  &  Latcha  C53,  Bolgov  &  Nikiforov  T63  and 
Heitkamper  C123  are  relevant. 


16 


In  the  acoustic  far-field  it  is  much  better  to  obtain  the 
transient  pressure  by  evaluating  numerically  a  Fourier  integral, 
obtained  via  a  stationary  phase  approximation,  by  an  application 
of  the  fast  Fourier  transform  algorithm.  The  Fortran  program 
which  computes  the  far-field  transient  radiation  is  small,  and 
does  not  include  any  proprietary  software.  Both  Sophie  Germain 
and  Timoshenko-Mindlin  plate  theories  are  included  in  the  program 
as  options,  and  'exact'  linear  theory,  of  layered  media,  using 
the  methods  of  Spicer  C9D,  could  be  added  without  too  much  dif¬ 
ficulty.  Thus,  a  comparison  of  predictions  from  'thin',  'thick' 
and  'exact'  theories  could  be  done  for  transient  sound  radiation, 
in  much  the  same  way  that  Clement  C133  has  done  so  for  the  time- 
harmonic  case. 

An  ad  hoc  treatment  of  the  contour  of  integration  used  to 
solve  the  Fourier  integral  has  shown  that  the  poles  in  the  integ¬ 
rand  have  a  physical  meaning,  viz.  direct  radiation  from  the 
excitation  point,  arriving  at  the  measuring  position  at  time 
t=R/c,  and  indirect  radiation  from  supersonic  waves,  coming  off 
the  surface  of  the  plate  at  oblique  angles,  arriving  before  t=R/c. 
These  features  merit  rigorous  investigation.  The  numerical 
results  shown  are  in  accordance  with  the  physical  interpretation. 
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FIG.  1  SECTION  OF  INFINITE  PLATE 
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FIG.  2  WARHAM.  EXCITATION  IS  A  UNIT  IMPULSE 
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FIG.  3  WARHAM.  UNIT  SIN  WAVE  SWITCHED  ON  FOR  5  CYCLES 

R«1.0m  9*45°  (a)  1kHz  (b)  10kHz 
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PRESSURE  IN  PASCALS 


FIG.  4  WARHAM.  UNIT  RECTANGULAR  PULSE  EXCITATION 

6«45°  (a)  R*0. 2m  (b)  R*1.0m. 
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FIG.  6  FFT.  EXCITATION  IS  A  UNIT  IMPULSE 

(a)  6*0°  <b)  0*45° 
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FIG.  7  FFT.  UNIT  SIN  WAVE  SWITCHED  ON  FOR  5  CYCLES 

(a)  1kHz  (b>  10kHz 
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FIG.  10  FFT.  COMPARISON  OF  PLATE  THEORIES 

UNIT  RECTANGULAR  PULSE  EXCITATION 
0=45°  (a)  THICK  <b)  THIN 
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FIG.  11  FFT.  COMPARISON  OF  PLATE  THEORIES 

UNIT  TRIANGULAR  PULSE  EXCITATION 

6=45°  (a)  THICK  (b)  THIN 
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